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Background: Pneumonia poses a significant risk in patients with moderate to severe chronic 
obstructive pulmonary disease but data are limited on the disease phenotypes most susceptible 
to pneumonia. 

Methods: Cluster analysis using a data-driven recursive partitioning algorithm was employed 
using baseline data from two pooled one-year randomized exacerbation trials (n=3,255) of 
fluticasone furoate/vilanterol or vilanterol alone to identify distinct patient groups at greatest 
risk of pneumonia or serious (hospitalization or death) pneumonia. 

Results: Five clusters were identified. Patients at greater risk of first pneumonia had more severe 
obstruction (forced expiratory volume in one second/forced vital capacity <46%) and either a 
body mass index <19 kg/m^ (hazard ratio 7.8, 95% confidence interval 4.7-13.0; n=144) or a 
pneumonia history and greater comorbidities (hazard ratio 4.8, 95% confidence interval 3.0-7.7; 
n=374) relative to the cluster with the lowest pneumonia risk (reference; n=1310). Multiple 
comorbidities and use of psychoanaleptics also contributed to an increased risk of pneumonia in 
more obstructed patients. Independent of cluster, use of inhaled corticosteroids was associated 
with pneumonia (hazard ratio 1 .89, 95% confidence interval 1 .25-2.84) and serious pneumonia 
(hazard ratio 2.92, 95% confidence interval 1.40-6.01). 

Conclusion: Cluster analysis can identify patient populations at risk for serious safety outcomes 
and inform risk management strategies to optimize patient management. The greatest risk for 
pneumonia was in subjects with multiple pneumonia risk factors. 

Keywords: chronic obstructive pulmonary disease, inhaled corticosteroids, long-acting 
p2-agonists, pneumonia, cluster analysis 

Introduction 

Community-acquired pneumonia results in greater morbidity and mortality in patients 
with chronic obstructive pulmonary disease (COPD) than in those without the disease.' 
The armual incidence of community-acquired pneumonia in the COPD population 
is estimated to be 22.4 episodes per 1,000 patients,^ and age, severity of COPD, and 
presence of comorbidity further increase the risk of hospitalization with community- 
acquired pneumonia. Findings further suggest an increase in pneumonia when COPD 
is treated with inhaled corticosteroids (ICS), although ICS also reduce the rates of 
moderate or severe exacerbations of COPD.' The mechanism by which ICS increase 
the risk of pneumonia is unclear, but may relate to reducing the inflammatory response, 
including neutrophils, in the respiratory tract.'' 

Fluticasone furoate/vilanterol (FF/VI) is a new once-daily ICS/long-acting 
P2-agonist (LABA) combination (FF/VI 100/25 |ig, BREO™ ELLIPTAtm, 
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GlaxoSmithKline, Research Triangle Park, NC, USA) that 
reduces exacerbations of COPD, improves lung function, 
and may increase the risk of pneumonia compared with the 
LABA (vilanterol) alone. ^ Factors indicated to be associated 
with this increased risk among patients with COPD treated 
with FF/VI are consistent with previous observations in the 
ICS/LABA class, and include older age, lower body mass 
index (BMI), current smoking, the occurrence of a previous 
pneumonia, and poorer lung function.^ '' 

We conducted a cluster analysis using data from two one- 
year studies of FF/VI and vilanterol to identify combinations 
of patient characteristics that place patients at greater risk of 
pneumonia for the purposes of risk management. 

Materials and methods 

Clinical study design and subjects 

The study designs, recruitment criteria, and procedures of the 
two studies analyzed herein have been previously reported.^ 
Briefly, patients aged &40 years with a documented diagnosis 
of COPD' and at least one moderate or severe exacerbation 
in the prior 12 months were randomized to 52 weeks of treat- 
ment with vilanterol 25 |J,g or FF/VI at strengths of 50/25, 
100/25, or 200/25 |J,g. The primary efficacy endpoint was the 
annual rate of moderate or severe exacerbations (requiring 
antibiotics and/or oral corticosteroids and/or resulting in 
hospitalization) in each treatment arm. Pneumonia was 
also reported and recorded. The diagnosis was based on the 
investigator's judgment based on the available clinical evi- 
dence, including a chest X-ray. Pneumonia was defined as 
an adverse event or as a serious adverse event if it resulted 
in hospitalization. 

Cluster analysis methodology 

Cluster analysis was employed using the R-package and 
RPART procedure r http://www.r-proiect.org/ ).^'^ Using 
the data-driven recursive partitioning algorithm to model 
time to first or severe pneumonia, the decision tree created 
clusters of patients based on the intent-to-treat population 
such that differences in risk of pneumonia were maximized 
among clusters. The size of the decision tree was deter- 
mined by minimizing the cross validation error; and cluster 
membership was then assigned to each patient based on 
the selected tree. The three doses of FF/VI were collapsed 
into one group to examine treatment versus nontreatment 
with FF/VI. Risk of pneumonia and serious pneumonia 
was calculated using Cox proportional hazards, with an 
initial model adjusted for study, smoking status, and region 



as a random effect. Because these variables were part of 
the initial model used to calculate risk, they could not be 
included as variables for defining the clusters. The results 
of the initial analysis did not identify treatment as a variable 
associated with clusters (see Results section), thus treat- 
ment was included as an explanatory variable within the 
final model. To avoid identification of spurious clusters, the 
number of patients within any one cluster was set at > 100. 
To avoid confounding by similar variables (eg, membership 
of a specific age group such as 65-75 years versus age as 
a continuous variable) each variable was assessed for col- 
linearity with every other variable using Pearson's correla- 
tion coefficients. Where two similar variables exhibited a 
Pearson's correlation coefficient value ^0.7 (eg, eosinophil 
count versus percentage) only one variable was retained 
based on the clinician's selection. Variables included in the 
analysis are noted in Table 1 . 

Hazard ratios (HR) and 95% confidence intervals (CI) 
were calculated for each cluster, using the cluster with 
the lowest incidence rate as the reference value. Baseline 
variables within each cluster are presented as proportions 
for categorical variables and medians (interquartile range) 
for continuous variables. Differences between variables in 
clusters were assessed using for categorical variables and 
the Wilcoxon rank sum test for continuous variables. 

Results 

Primary efficacy and safety 

The primary data of these studies have already been reported.^ 
More pneumonia, serious pneumonia, and fatal pneumonia 
events occurred with FF/VI than vilanterol (Table 2). 

Cluster analysis 

This analysis yielded five clusters combining clinical and 
demographic attributes, including forced expiratory volume 
in one second/forced vital capacity (FEV^FVC) ratio, psy- 
choanaleptic use, BMI, and history of pneumonia (Figure 1 ). 
The final tree was based on combining the similar trees 
determined for all pneumonia events (Figure S 1 ) and serious 
pneumonia events (subset of pneumonia events resulting in 
hospitalization; Figure S2). Treatment (FF/VI or vilanterol) 
was not identified as a variable in the tree that maximized 
the difference in risk between clusters, so was included in 
the final Cox model. The incidence of pneumonia and seri- 
ous pneumonia was lowest in cluster 1 (2.52% and 0.92%, 
respectively) during the one-year study period and highest 
in cluster 5 (18.75% and 9.03%, respectively), generally 
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Table I Characteristics and variables considered in initial cluster 
analysis 



Characteristic/variable: class Characteristic/variable: specific 



Spirometry 


• Post-bronchodilator FEV| L 




• Post-bronchodilator FEV|% predicted 




• FEV|% reversibility 




• Post-bronchodilator FEV|/FVC ratio 




• Post-bronchodilator FVC% predicted 


Demographics 


• Age (years) 




• Sex 




• BMI (kg/m^) 




• Smoking (current/former) 




• Smoking pack-years 


Race 


• African American/African heritage 




• American Indian/Alaska Native 




• Asian 




• Mixed 




• White 


COPD history 


• Duration of COPD (years) 




• COPD type (chronic bronchitis 




and/or emphysema) 


History of exacerbation 


• < 1 versus £2 requiring oral 


(past year) 


corticosteroids and/or antibiotics 




• < 1 versus a2 requiring 




hospitalization 0, 1,2, >2 or more 


Medical history 


• Targeted medical conditions* 




• Vaccination for influenza 




(yes versus no) 




• Vaccination for pneumonia 




(yes versus no) 


Medications at baseline/initiated 


• ATC classification** 


prior to randomization 




Laboratory investigations 


• Blood urea nitrogen 




• Eosinophils 




• Hemoglobin 




• Lymphocytes 




• Neutrophils 




• Red blood cells 




• White blood cells 


Randomized treatment 


. VI 25 ng 




• FF/VI treatment (50/25 \xg, 




100/25 |ig, 200/25 |ig) 



Notes: ^Targeted medical conditions: cardiovascular history/risk, vascular disorders, 
metabolism disorders, history of pneumonia, cardiac disorders, eye disorders, arrhythmia; 
^^medications were coded based on the ATC of interest. Some medications have more 
than one indication and were categorized in more than one therapeutic class (eg, aspirin 
for pain relief versus cardiovascular disease prevention). Because the reasons for taking 
the medications were not collected systematically, the classification was based on the 
ATC classes and not on the underlying indication for an individual patient. It is possible 
that some medications were counted in more than one ATC category. Medication classes 
included in the analysis included: antihistamines, antihypertensives, antithrombotics, 
diuretics, diabetes medications, psychoanaleptics, psycholeptics, agents acting on the 
renin-angiotensin system {eg, angiotensin-converting enzyme inhibitors, angiotensin II 
agonists), antihemorrhagics, antianemic preparations, beta-blockers, cardiac therapies, 
calcium channel blockers, lipid-modifying agents, anti-inflammatory and antirheumatic 
products, medications for bone diseases (including muscle pain), antihemorrhagics, and 
vasodilators. 

Abbreviations: ATC, Anatomic Therapeutic Category; COPD, chronic obstructive 
pulmonary disease; BMI, body mass index; FEV|, forced expiratory volume in 
one second; FVC, forced vital capacity; FF/VI, fluticasone furoate/vilanterol; 
VI, vilanterol. 



Table 2 Pneumonia, 


serious 


pneumonia, 


and fatal 


pneumonia 


events with FF/VI and VI 








n (%) 


FF/VI 


r r/ VI 


rr/ VI 


VI 




50 ng/ 


1 nn no/ 




IICT 




25 Hg 


" Jig 




/rt— A 1 n\ 
^^n— o 1 o f 




(n=820) 


in — ouo 1 


(n=8 1 1 ) 




Pneumonia 










Subjects 


48 (5.9) 


51 (6.3) 


55 (6.8) 


27 (3.3) 


Events 


54 


58 


65 


28 


Serious pneumonia 










Subjects and events 


24 (2.9) 


25 (3.1) 


23 (2.8) 


8(<l) 


Fatal pneumonia 










Subjects and events 


0 


1 


7* 


0 



Note: *One subject in the 200/25 treatment group experienced a fatal exacerbation 
of chronic obstructive pulmonary disease, and was noted to have concurrent 
pneumonia at the time of death. 

Abbreviations: FF/VI, fluticasone furoate/vilanterol: VI, vilanterol. 



increasing from left to right (Figure 2A and C). Relative to 
cluster 1 (reference), HRs for the time to first pneumonia in 
clusters 2-5 ranged from a 2.2-fold increased risk of first 
pneumonia in cluster 3 versus cluster 1 to a 7.8-fold increased 
risk of first serious pneumonia in cluster 5 versus cluster 1 
(Figure 2B and D). 

Patient clusters 

The clusters identified comprised the following: patients 
with poorer lung function (FEVj/FVC <46.05%) and lower 
BMI < 19 kg/m^ (cluster 5); patients with poorer lung func- 
tion, normal BMI and prior history (cluster 4) or no history 
(cluster 3) of pneumonia; and patients with better lung 
function in receipt (cluster 2) or not in receipt (cluster 1) 
of psychoanaleptics (89% antidepressants). There were 
significant baseline differences between clusters assessed 
by demographics, spirometry, COPD history, or COPD type 
(Table 3) and comorbidities or concomitant medications 
(Table 4). There were no clinically significant differences 
in laboratory measures across clusters (data not shown). 

Patients with the greatest risk of pneumonia and seri- 
ous pneumonia (cluster 5) predominantly had a history of 
emphysema, the highest prevalence of current smoking, 
the greatest extent of reversibility, and were more likely to 
have experienced at least two exacerbations in the prior year 
compared with those in other clusters. Patients in cluster 5 
also had considerably fewer comorbidities than other clusters 
and were in receipt of the fewest concomitant medications. 
Patients in clusters 4 and 2 showed higher rates of comorbid- 
ity and concomitant medication. 

Patients in cluster 3 were predominantly male and had 
COPD for the greatest amount of time, but were the least likely 
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FEV/FVC ratio (%) 



>46.05 



<46.05 



Psychoanaleptic use 



No 



BMI (kg/m2) 



>19 



Yes 



<19 



Pneumonia history 



No 



Yes 



Incidence rate, % 



^^i^^^ ^^m^^ ^^^^^^ ^^m^^ ^^^^^^ 

^^^^^^ ^^^^^^ ^gm^^ ^^^^^^ ^QQi^ 



1st 

pneumonia 


FF/VI 


2.8 


9.3 


5.7 


12.5 


21.3 


VI 


1.6 


3.6 


3.5 


5.8 


11.1 


1st serious 
pneumonia 


FF/VI 


1.0 


3.9 


2.4 


8.0 


10.2 


VI 


0.6 


0 


1.0 


1.2 


5.6 



Figure I Cluster analysis tree of first and first serious pneumonia. 

Abbreviations: BMI, bone mass index; FEV|, forced expiratory volume in one second; FVC, forced vital capacity; FF/VI, fluticasone furoate/vilanterol; VI, vilanterol. 
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Figure 2 Incidence (A and C) and hazard ratios {B and D) for first pneumonia (A and B) and first serious pneumonia (C and D) by cluster. 
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Table 3 Patient characteristics by cluster 



Parameter 



Cluster I 
(less 

obstruction) 
n=l,3IO 



Cluster 2 

(less obstruction, 

psychoanaleptic 

users, high 

comorbidity) 

n=238 



Cluster 3 

(greater 

obstruction, 

no pneumonia 

history) 

n=l,l89 



Cluster 4 
(greater 
obstruction, 
high comorbidity, 
pneumonia history) 
n=374 



Cluster 5 
(greater 
obstruction, 
low BMI) 
n=l44 



P-value 



Demographics 

Age, median (IQR) 

Male, n (%) 

Body mass index at 

screening, median (IQR) 
Ethnicity, n (%) 

Not Hispanic or Latino 

Hispanic or Latino 
Smoking status 

Former, n (%) 

Current, n (%) 

Number of pack years, 

median (IQR) 
Spirometry 

% predicted post- 

bronchodilator FEV| 

at baseline, median (IQR) 

Percent reversibility (%) 

at screening, median (IQR) 

Reversible*, n (%) 
Duration of COPD, n (%) 

<S years 

5 to < 1 0 years 

1 0 to < 1 5 years 

a I S years 



64 (56-70) 
701 (53.5) 
27.3 (24.1-31.3) 



1,000 (76.3) 
310 (23.7) 

736 (56.2) 
574 (43.8) 
38.0 (22.0-50.0) 



60 (55-67) 
67 (28.2) 

29.6 (25.0-34.1) 



224 (94.1) 

14 (5.9) 

99 (41.6) 
139 (58.4) 
46.0 (35.0-60.0) 



65 (59-71) 
788 (66.3) 

25.6 (23.0-28.6) 



986 (82.9) 
203 (17.1) 

702 (59.0) 
487 (41.0) 
42.0 (30.0-60.0) 



65 (59-71) 

234 (62.6) 
25.7 (22.6-29.0) 



328 (87.7) 
46 (12.3) 

224 (59.9) 

150 (40.1) 
45.0 (32.0-64.0) 



5 1 .9 (43.9-60.3) 52. 1 (43.7-60.2) 35.0 (27.2-43.7) 33.5 (26.0-4 1 .6) 



9.0 (2.4-17.5) 

372 (28.5) 

581 (44.4) 
432 (33.0) 
158 (12.1) 
139 (10.6) 



1 3.6 (6.9-22.6) 
101 (42.6) 

I 1 3 (47.5) 

65 (27.3) 
38 (16.0) 
22 (9.2) 



Exacerbation requiring OCS/antibiotic in past 12 months, n (%) 

0 94 (7.2) 16 (6.7) 

1 865 (66.0) 140 (58.8) 

2 248(18.9) 50(21.0) 
>3 103 (7.9) 32 (13.4) 

Exacerbation requiring hospitalization in past 1 2 months, n (%) 



0 
I 

>2 

COPD type, n (%) 
History of bronchitis 
History of emphysema 



1,031 (78.7) 
219 (16.7) 
60 (4.6) 

947 (72.8) 

582 (44.8) 



195 (81.9) 
36 (15.1) 
7 (2.9) 

170 (71.4) 
1 28 (53.8) 



14.3 (6.0-25.2) 

362 (3 1 .2) 

466 (39.2) 
347 (29.2) 
204 (17.2) 
172 (14.5) 

81 (6.8) 
722 (60.7) 
279 (23.5) 
107 (9.0) 

1,001 (84.2) 

157 (13.2) 
31 (2.6) 

756 (63.9) 
729 (61.6) 



16.2 (8.6-25.4) 

101 (27.6) 

I 17 (31.3) 
128 (34.2) 
82 (21.9) 
47 (12.6) 

40 (10.7) 
196 (52.4) 
82 (21.9) 
56 (15.0) 

256 (68.4) 
104 (27.8) 
14(3.7) 

214 (57.2) 
260 (69.5) 



64(56-70) <0.000l 
80 (55.6) <0.000l 
17.6(16.4-18.3) <0.000l 



1 33 (92.4) 
I I (7.6) 



36 (25.2) 

56 (38.9) 
56 (38.9) 

17 (I 1.8) 
15 (10.4) 

19 (13.2) 

69 (47.9) 
34 (23.6) 
22 (15.3) 

I 10 (76.4) 
31 (21.5) 
3(2.1) 

69 (47.9) 
I 10 (76.4) 



<0.000l 



<0.000l 



55 (38.2) 
89(61.8) 

43.0 (30.0-56.5) <0.000l 



28.0 (21.8-34.9) <0.000l 



18.1 (8.3-28.8) <0.000l 



0.0002 
<0.000l 



<0.000l 



<0.000l 



<0.000l 
<0.000l 



Notes: *ReversibiIity indicated by ^ 12% and 200 mL change in lung function following bronchodilator. The most prevalent/greatest parameter across clusters is shown in bold. 
Abbreviations: BMI, bone mass index; COPD, chronic obstructive pulmonary disease; OCS, oral corticosteroids; IQR, interquartile range; FEV|, forced expiratory volume 
in one second. 



to have experienced an exacerbation requiring hospitalization 
in the prior year; this cluster was primarily differentiated 
from the other clusters on the basis of poorer lung function 
and the absence of other defining characteristics. 

Cluster 2 subjects were predominantly female, in addition 
to having a higher prevalence of comorbidity (and typically 
the highest receipt of concomitant medication). This cluster 
contained the fewest Hispanic and Latino patients, had the 
least impaired lung function, greatest proportion of current 



smokers, and greatest proportion of reversible patients. 
Forty-three percent of subjects in cluster 2 were also using 
psycholeptic medications, over half of whom were exposed 
to benzodiazepines (Table S 1 ). 

Cluster 1 (the reference population) contained the great- 
est number of patients overall, and were more likely to be 
Hispanic or Latino and more likely to have a history of 
COPD characterized by bronchitis than patients in the other 
clusters. 
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Table 4 Patient comorbidities and concomitant medications by cluster 



Parameter, n (%) 


Cluster 1 


Cluster 2 


Cluster 3 


Cluster 4 


Cluster 5 


P-value 




(less 


(less obstruction, 


(greater 


(greater 


(greater 






obstruction) 


psychoanaleptic 


obstruction. 


obstruction. 


obstruction. 






n=l,3 10 


users, high 


no pneumonia 


high comorbidity. 


low BMI) 








connorbidity) 


history) 


pneumonia history) 


n— 1 '•'1 








n=238 


n=l , 1 89 


n=374 






Comorbidities 














Cardiovascular history/risk 


806 (6 1 .5) 


1 79 (75.2) 


706 (59.4) 


253 (67.6) 


65 (45. 1 ) 


<0.000l 


Vascular disorders 


63 1 (48.2) 


1 34 (56.3) 


540 (45.4) 


1 96 (52.4) 


40 (27.8) 


<0.000l 


Metabolism disorders 


505 (38.5) 


143 (60.1) 


413 (34.7) 


193 (51.6) 


40 (27.8) 


<0.000l 


History of pneumonia 


241 (18.4) 


87 (36.6) 


0 (0.0) 


374 (100) 


38 (26.4) 


<0.000l 


Cardiac disorders 


220 (16.8) 


48 (20.2) 


162 (13.6) 


74 (19.8) 


22 (15.3) 


0.0136 


Eye disorders 


140 (10.7) 


31 (13.0) 


130 (10.9) 


84 (22.5) 


22 (15.3) 


<0.000l 


Arrhythmia 


80 (6.1) 


13(5.5) 


38 (3.2) 


21 (5.6) 


1 0 (6.9) 


0.0 1 07 


Concomitant medications 














On-treatment pneumonia 


1 1 7 (8.9) 


15(6.3) 


51 (4.3) 


38 (10.2) 


6 (4.2) 


<0.000l 


Vaccine 














On-treatment influenza vaccine 


312 (23.8) 


55 (23.1) 


220 (18.5) 


1 19 (31.8) 


27 (18.8) 


<0.000l 


Agents acting on the 


419 (32.0) 


85 (35.7) 


365 (30.7) 


122 (32.6) 


18 (12.5) 


<0.000l 


renin-angiotensin system 














Anti-inflammatory and 


406 (3 1 .0) 


127 (53.4) 


388 (32.6) 


144 (38.5) 


29 (20.1) 


<0.000l 


antirheumatic products 














Lipid-modifying agents 


342 (26. 1 ) 


1 10 (46.2) 


279 (23.5) 


1 28 (34.2) 


22 (15.3) 


<0.000l 


Antithrombotics 


308 (23.5) 


78 (32.8) 


308 (25.9) 


1 1 2 (29.9) 


25 (17.4) 


0.00 1 1 


Diuretics 


256 (19.5) 


69 (29.0) 


238 (20.0) 


79 (2 1 . 1 ) 


19 (13.2) 


0.0027 


Beta-blockers 


151 (1 1.5) 


33 (13.9) 


88 (7.4) 


35 (9.4) 


5 (3.5) 


0.0001 


Diabetes medications 


1 46 ( 1 1 . 1 ) 


38 (16.0) 


92 (7.7) 


39 (10.4) 


1 (0.7) 


<0.000l 


Psycholeptics 


123 (9.4) 


102 (42.9) 


142 (1 1.9) 


65 (17.4) 


18 (12.5) 


<0.000l 


Antihistamines 


1 16 (8.9) 


53 (22.3) 


1 1 1 (9.3) 


53 (14.2) 


10 (6.9) 


<0.000l 


Antihypertensives 


53 (4.0) 


6 (2.5) 


51 (4.3) 


31 (8.3) 


8 (5.6) 


0.0038 


Antihemorrhagics/ 


48 (3.7) 


30 (12.6) 


43 (3.6) 


21 (5.6) 


4 (2.8) 


<0.000l 


antianemia preparations 














Drugs for bone diseases 


39 (3.0) 


15(6.3) 


37(3.1) 


29 (7.8) 


1 1 (7.6) 


<0.000l 


(including muscle pain) 














Psychoanaleptics 


0 (0.0) 


238 (100) 


172 (14.5) 


66 (17.6) 


23 (16.0) 


<0.000l 


Note: Most prevalent parameter across 


clusters is shown in 


bold. 











Abbreviation: BMI, body mass index. 

Treatment effects 

In each cluster, pneumonia or serious pneumonia occurred 
more frequently with any strength of FFA'^I therapy than with 
vilanterol alone (see Table in Figure 1). Across all clusters 
and relative to treatment with vilanterol alone, treatment with 
FF/VI at any strength (ie, all strengths combined) was associ- 
ated with a significantly greater risk of first pneumonia (HR 
1.89, CI 1.25-2.84) and first serious pneumonia (HR 2.92, 
CI 1.40-6.01). No interaction was found between treatment 
and cluster for first pneumonia (/'=0.9741) or first serious 
pneumonia (/'=0.8089). Doses of FF/VI were balanced across 
the clusters (/'=0.2064, Table S2). 

Effects of variables used in the model 

Smoking status (current) as a main effect was not associated 
with a significantly increased risk of first pneumonia (HR 1 .0, 



CI 0.7-1.4) or first serious pneumonia (HR 0.9, CI 0.6-1.5). 
The distribution of geographic region across clusters was sig- 
nificantly different (P^O.OOO 1 ) relative to the total distribution 
(Table 5). Overall, approximately one third of patients were 
from the USA, a quarter were from Europe, and the remainder 
were balanced between the two "other" regions (Argentina, 
Chile, Mexico, Peru, and the Philippines as other region 1, 
and Australia, Canada, and South Africa as other region 2). 
In cluster 2 (in receipt of psychoanaleptics, higher lung fiinc- 
tion, more comorbidity, predominantly female), 61.3% of the 
population were from the USA, while compared with the total 
population proportionately fewer patients were from Europe 
or "other" region 1. The USA also provided the majority of 
patients in cluster 4 (poorer lung function, higher proportions 
of comorbidity, prior pneumonia) and cluster 5 (poorer lung 
function, lower BMI, greatest pneumonia risk). 
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Table 5 Geographic contribution to total population and each cluster 



Region, % 


Total 


Cluster 1 


Cluster 2 


Cluster 3 


Cluster 4 


Cluster 5 


patients 


n=3,255 


(less 


(less obstruction, 


(greater obstruction. 


(greater obstruction. 


(greater 






obstruction) 


psychoanaleptic users. 


no pneumonia 


high comorbidity. 


obstruction. 






n=l ,3 1 0 


high comorbidity) 


history) 


pneumonia history) 


low BMI) 








n=238 


n=l,l89 


n=374 


n=l44 


USA 


34.0 


26.6 


61.3 


30.3 


53.7 


36.8 


Europe 


26.1 


27.7 


14.3 


30.5 


18.2 


13.9 


Other 1 


21.4 


28.8 


6.3 


19.2 


1 1.8 


22.2 


Other 2 


18.5 


17.0 


18.1 


20.0 


16.3 


27.1 



Notes: Other I, Argentina, Chile, Mexico, Peru, Philippines; Other 2, Australia, Canada, South Africa. The most prevalent region in each cluster is shown in bold. 
Abbreviation: BMI, body mass index. 



Discussion 

Cluster analysis of two replicate one-year studies of FF/VI 
versus vilanterol identified five clusters of patients with 
differing risks for occurrence of pneumonia or serious 
pneumonia. Using the cluster with the lowest incidence of 
pneumonia (cluster 1) as a reference, the risk of pneumonia 
or serious pneumonia was significantly increased in each of 
the remaining clusters. The clusters identified were defined 
by combinations of the established pneumonia risk factors 
of lung function, BMI, and prior history of pneumonia, and 
also by the use of psychoanaleptic medication. FEVj/FVC 
was the main characteristic used to describe the clusters, and 
clusters 3, 4, and 5 were characterized by greater obstruction, 
as defined by FEV/FVC, than clusters 1 and 2. 

Cluster 5 had the highest risk of pneumonia, and subjects in 
this cluster showed the greatest extent of lung fiinction impair- 
ment, defined as FEVj% predicted, and the lowest BMI of all 
the clusters. Cluster 4 subjects had more impaired lung function 
(FEVj% predicted), more comorbid conditions, and a prior 
history of pneumonia compared with cluster 1, while those 
in cluster 3 had only more impaired lung fianction (FEVj% 
predicted) compared with cluster 1. Cluster 2 subjects also 
showed an increased risk of pneumonia relative to cluster 1 ; 
cluster 2 subjects differed from those in cluster 1 by use of 
psychoanaleptics and an increased prevalence of comorbidi- 
ties. The algorithm did not identify treatment with FFA^I versus 
vilanterol as a variable maximizing the difference in risk of 
pneumonia between clusters; however, treatment with FFA^^I 
versus vilanterol alone did result in a significantly increased 
risk of pneumonia or serious pneumonia in each cluster There 
was no differential treatment effect across clusters. 

Vanfletern et aP" have identified different COPD pheno- 
types based on the number and types of comorbidities, and 
the phenotypes identified loosely track with the clusters of 
varying pneumonia risk identified in our analysis, includ- 
ing low comorbidity (cluster 1), psychological (cluster 2), 



cardiovascular (clusters 2 and 4), metabolic (clusters 2 
and 4), and cachectic (cluster 5) groups. In our analysis, the 
initial split of the population occurred based on lung fiinction 
(FEVj/FVC ratio <46%). Poorer lung function is known to 
increase the risk of pneumonia in patients with COPD,' " 
so it is perhaps unsurprising that our analysis defined lung 
function as the initial variable associated with an increased/ 
decreased risk of pneumonia. Increasing obstruction mea- 
sured as FEVj and increasing prevalence of emphysematous 
COPD generally tracked with increasing risk among the 
cohorts. 

The poorer lung fiinction population subsequently split on 
the basis of BMI and patients with both poorer lung fiinction 
and a BMI <19 kg/m^ formed the cluster (cluster 5) with 
the greatest risk and incidence of pneumonia and serious 
pneumonia. This finding in part confirms that of the initial 
analysis, in which BMI was associated with the occurrence of 
pneumonia,'^ and also aligns with other findings." It follows 
that in patients with both poor lung function and low BMI, 
the risk of pneumonia would be increased as a consequence 
of the presence of both risk factors. 

Lower lung function, a prior pneumonia event, and 
greater prevalence of comorbidities were the defining charac- 
teristics of cluster 4. In COPD, it is known that the incidence 
of community-acquired pneumonia is approximately doubled 
in patients with a prior event,^ which suggests that, as has 
been shown for exacerbations of COPD," the occurrence of 
a prior event may increase the risk of a recurrent event. 

The finding that cluster 2 comprised patients with better 
lung function (FEVj/FVC >46%), more psychoanaleptic 
use, and frequent psycholeptic use suggests that receipt 
of psychotropic medications presents an increased risk of 
pneumonia in COPD that is distinct from established fac- 
tors. Other studies have reported that both psychotropic 
drugs"* and the underlying psychiatric/neurological disorders 
they treat'^"' increase the risk of pneumonia. Patients with 
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psychiatric and neurological disorders can have impairment 
in both the perception of symptoms and the seeking of health 
care, and may be more likely to smoke." Patients with COPD 
who smoke have a higher prevalence of depression and are 
at greater risk of developing depression than those who do 
not.'* Depressive symptoms are known to pose a greater risk 
of exacerbation, exacerbation-related hospitalization, and 
death in COPD." 

Patients taking psychotropic medications are also at 
risk for aspiration pneumonia because of dysphagia"* '^ 
and relaxation of lower esophageal sphincter tone that can 
aggravate gastroesophageal reflux.^" The most common class 
of psychotropic drugs used by patients in these studies was 
benzodiazepines, which were estimated to be dispensed to 
nearly one third of patients with COPD in a Canadian study.^' 
A recent nested case-control and survival analysis in a popu- 
lation-based cohort indicated that benzodiazepines increase 
the risk of pneumonia by nearly 50%,'' although an earlier 
study suggested that this risk was related to concomitant 
opioid use.^^ Aside from the aspiration risk cited, it has been 
suggested that benzodiazepine-related pneumonia may be 
due to immunosuppression that is mediated via activation of 
gamma-aminobutyric acid A receptors on macrophages.^" 

Another important aspect related to increased risk 
of pneumonia identified in this analysis is the presence 
of multiple comorbidities. The presence of a history or risk of 
cardiovascular disease and the incidence of coronary artery 
disease and cardiac, metabolic, or vascular disorders was 
highest in clusters 2 and 4, in whom the risk of pneumonia 
was greater than in both cluster 1 (the reference population) 
and cluster 3. Both congestive heart failure and peripheral 
vascular disease are independent risk factors for the occur- 
rence of community-acquired pneumonia.^ Similarly, both 
cardiovascular disease and diabetes place patients with 
COPD at an increased risk of hospitalization for pneumonia. 
This would suggest that an increased risk of pneumonia 
would be most likely to be present in subjects with comorbid 
cardiovascular disease treated with FF/VI versus vilanterol. 
However, in other analyses of the data utilized here (not 
shown), subjects with comorbid cardiovascular disease 
treated with FF/VI had a numerically lower risk for occur- 
rence of pneumonia when compared with subjects with no 
comorbid cardiovascular disease treated with vilanterol. The 
definition of comorbid cardiovascular disease used in this 
analysis was a prior or current diagnosis of coronary artery 
disease, myocardial infarction, arrhythmia, congestive heart 
failure, hypertension, cerebrovascular accident, diabetes 
mellitus, or hypercholesterolemia. This broad definition 



may have had a confounding effect by classing together 
morbidities, and perhaps does not adequately reflect the 
complexity of COPD (eg, a single factor versus multiple 
factors contributing to risk of pneumonia). 

It is also notable that only a minority of patients were 
vaccinated for pneumonia (range across clusters, 4%-9%) or 
influenza (range 19%-23%). Vaccination rates were highest 
for cluster 4, which had a higher level of obstruction, pneu- 
monia history, and more frequent comorbidities relative to 
other clusters. Cluster 4 patients included disproportionate 
numbers from the USA, suggesting that greater vaccination 
rates and comorbidities in this cluster could relate to differ- 
ences in health systems and lifestyle or nutrition relative to 
other regions. Imbalances in geographic distribution across 
other clusters support this possibility. The low vaccination 
rates in these COPD study patients suggests a modifiable 
risk of pneumonia exists that can be further addressed within 
clinical practice settings and by vaccine uptake strategies. 
Although the evidence is not consistent, the GOLD (Global 
Initiative for Chronic Obstructive Lung Disease) consensus 
document^^ and Centers for Disease Control and Prevention^* 
recommend that influenza and pneumonia vaccination be 
offered to every COPD patient, where vaccination is more 
effective in older patients, severe COPD patients, and those 
with comorbid cardiovascular disease. Despite the increased 
risk of pneumonia in these groups, the patients with greatest 
risk of pneumonia are also those who receive the greatest 
benefit of ICS, including those with increased obstruction 
and low BMI.'' 

Treatment allocation was balanced across clusters, which 
would be expected due to randomization. There was a small 
imbalance in cluster 5, with a slight underrepresentation of 
FF/VI 100/25 iig (17.4%) and overrepresentation of FF/VI 
200/25 jig (31.3%), compared with the 25% expected. This 
variation likely resulted due to chance because this was the 
smallest of the clusters (n=144), but could have contributed 
to the observed increased risk of pneumonia. Based on the 
benefit-risk profile,^ the 200/25 dose of FF/VI was not pro- 
gressed for registration in COPD. 

Regardless of cluster, there is an approximate two-fold 
increase in the relative risk of pneumonia and a nearly three- 
fold increase of serious pneumonia among patients treated 
with FF/VI relative to vilanterol; however, the absolute risk 
of pneumonia is greatest among those who benefit most from 
FF/VI (eg, clusters 4 and 5).' Risk factors that put individu- 
als at greatest risk of pneumonia (eg, low BMI, increased 
airflow obstruction, multiple morbidities, and prior history 
of pneumonia) can be collectively assessed when making 
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treatment decisions. If the greatest risk of pneumonia had 
been found among patients with modest benefit from FF/VI, 
the benefit-risk profile may merit a bronchodilator alone 
(eg, vilanterol) rather than FF/VI, to maintain a positive 
benefit-risk profile. 

Cluster analysis is a usefiil tool to examine heterogeneous 
diseases such as COPD^' and has been employed to aid in the 
identification of distinct clusters or phenotypes of COPD,^^ 
including those at greater risk of mortality.^'-'" Interestingly, 
there is some overlap in the phenotypes identified as being 
at greater risk of mortality and pneumonia, including those 
with multiple comorbidities and moderate to severe airflow 
limitation^'-^" and those with severe airflow limitation with 
low BMI.^' Cluster analysis differs principally from sub- 
group analysis in that it seeks to first maximize differences 
between groups of patients using a data-driven algorithm 
and then identify multiple characteristics which define those 
patients. In comparison, subgroup analysis first defines 
patients by individual characteristics one at a time. 

This analysis benefited from a large dataset representing 
a broad international patient population, which is an impor- 
tant element in allowing robust clusters to be identified. 
However, the generalizability of these findings is limited 
to patients with a history of COPD exacerbation in the past 
year. An additional limitation is the lack of validation via a 
second dataset or split sample; however, results are consis- 
tent with what is known about risk factors for pneumonia. 
Given the small number of pneumonia events, we were 
not able to evaluate the effect of ICS dose on pneumonia 
within clusters. Finally, definitions of pneumonia employed 
in the original studies and in this analysis did not include 
any standardized assessment of pneumonia severity, such 
as CURB-65 or Pneumonia Severity Index.^' The definition 
of serious pneumonia was based on hospitalization due to 
pneumonia, which may vary across health systems, relative 
to standardized criteria. 

The present analysis identified five clusters of patients 
with differing risks of pneumonia occurrence over one year 
of therapy with FF/VI or vilanterol. The clusters were in 
part defined by well established risk factors, such as BMI 
and lung function and as such, the outcomes can be consid- 
ered confirmatory of other findings. The clusters were also 
defined by less established risk factors, such as prior occur- 
rence of pneumonia, receipt of psychotropic medications, 
and presence of multiple comorbidities, and as such can be 
considered hypothesis-generating. Our analysis suggests 
that physicians should consider many factors and combina- 
tion of risk factors when treating COPD patients with an 



ICS/LABA such as FF/VI, balancing the risks and benefit 
in the individual patient. 
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Supplementary materials 



Table SI Psychoanaleptic and psycholeptic prescription 


in the five clusters 








n ^/o OT Cluster 


^liic^ai- 1 

^lusner 1 


Cluster 2 


Cluster 3 


Cluster 4 


Cluster 5 


populations 




(less obstruction, 


(greater 


(greater 


(greater 




obstruction^ 


psychoanaleptic 


obstruction. 


obstruction. 


obstruction. 




n— 1 1 J 1 u 


users, high 


no pneumonia 


high comorbidity. 


low BMI) 






comorbidity) 


history) 


pneumonia history) 


n=l44 






n=238 


n=l,l89 


n=374 




Psychoanaleptics 












Any 


0 (0.0) 


zJo ( 1 00) 


172 (14.5) 


66 ( 1 7.6) 


23 (16.0) 


Antidepressants 


0 (0.0) 


1 1 T /OQ 1 \ 

2 1 2 (o7. 1 ) 


14/ (12.4) 


64 ( 1 7. 1 ) 


22 (15.3) 


Psychostimulants for ADHD/nootropics 


0 (0.0) 


ii (1 J.V) 


LI (2.J) 




1 ir\ T\ 
1 (0./) 


Antidementia medication 


0 (0.0) 


7 (2.9) 


4 (0.3) 


2 (0.5) 


0 (0.0) 


Psycholeptics 












Any 


1 23 (9.4) 


102 (42.9) 


142 (1 1.9) 


65 (17.4) 


18 (12.5) 


Antianxiety medication 


82 (6.3) 


68 (28.6) 


100 (8.4) 


38 (10.2) 


1 3 (9.0) 


Hypnotics/sedatives 


S3 (4.0) 


40 (16.8) 


56 (4.7) 


31 (8.3) 


6 (4.2) 


Antipsychotics 


10 (0.8) 


2! (8.8) 


12 (1.0) 


9 (2.4) 


1 (0.7) 


Benzodiazepine/derivative 


70 (5.3) 


59 (24.8) 


87 (7.3) 


34 (9.1) 


9 (6.3) 



Abbreviation: ADHD, attention-deficit hyperactivity disorder; BMI, body nnass index. 



FEV7FVC ratio (%1 



>46.05 



Psychoanaleptic use 



No 



<46.05 



BMI (kg/m^ 



>18.78 



Yes 



<18.78 



Pneumonia history 



No 



Yes 




Incidence rate, % 



2.52 



7.98 



5.18 



10.81 



19.85 



Figure SI Cluster analysis tree of first pneumonia. 

Abbreviations: BMI, bone mass index; FEV|, forced expiratory volume in one second. 
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Table S2 Distribution of treatment by cluster 



Treatment 


Total 


Cluster 1 


Cluster 2 


Cluster 3 


Cluster 4 


Cluster 5 


(% of patients) 


N=3255 


(less 


(less obstruction, 


(greater 


(greater 


(greater 






obstruction) 


psychoanaleptic users, obstruction. 


obstruction. 


obstruction. 






n=l3IO 


high comorbidity) 


no pneumonia 


high comorbidity. 


low BMI) 








n=238 


history) 


pneumonia history) 


n=l44 










n=l 189 


n=374 




VI 25 


25.1 


24.7 


23.5 


26.6 


23.3 


25.0 


FF/VI 50/25 


25.2 


24.7 


27.3 


23.5 


30.5 


26.4 


FF/VI 100/25 


24.8 


25.7 


25.6 


24.6 


24.3 


17.4 


FF/VI 200/25 


24.9 


25.0 


23.5 


25.3 


22.0 


31.3 



Abbreviations: BMI, body mass index; FF/VI, fluticasone furoate/vilanterol; VI, vilanterol. 



FEV % predicted 



>41.95 



>19.38 




<41.95 



BMI (kg/m^; 



<19.38 



Pneumon 


a history 


No 


Yes 






Incidence rate, % 



1.22 



1.96 



6.49 



9.47 



Figure S2 Cluster analysis tree of first serious pneumonia. 

Abbreviations: FEV|, forced expiratory volume in one second; BMI, body mass index. 
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